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Atomic clusters of various sizes irradiated with short intense pulses 
of VUV radiation 
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Abstract. - Non-equilibrium processes following the irradiation of atomic clusters with short 
pulses of vacuum ultraviolet radiation are modelled using kinetic Boltzmann equations. The 
dependence of the ionization dynamics on the cluster size is investigated. The predictions on: (i) 
the maximal and average ion charge created, (ii) ion charge state distribution, (iii) average energy 
absorbed per atom, (iv) spatial charge distribution, and (v) thermalization scales are obtained for 
spherical xenon clusters containing: 20, 70, 2500 and 90000 atoms. These clusters were exposed to 
single rectangular pulses of vacuum ultraviolet radiation of various pulse intensities, / ~ 10 12 — 10 14 
W/cm 2 and durations, At < 50 fs, at a fixed integrated radiation flux of F = / • At = 0.4 J/cm 2 . 
The results obtained are found to be in good agreement with the available experimental data, 
especially the dependence on the cluster size, if it is assumed that the ions from the positively 
charged outer layer of the cluster constitute the dominant contribution to the experimentally 
measured ion charge state distribution. 



Atomic clusters are physical objects of nanometer size, 
consisting of a few up to tens of thousands of closely 
packed atoms. Their physical properties put them onto 
the border between solid state and gas phase. Clusters are 
excellent objects to test the dynamics within samples irra- 
diated with short wavelength radiation from free-electron- 
lasers (FELs) [1-4]. By studying the size dependence the 
influence of the atomic and condensed matter effects on 
the ionization dynamics can be explored. Cluster studies 
are important for planned experiments with FELs within 
the solid state physics, material sciences and for studies 
of the extreme states of matter [5]. Accurate predictions 
on the ionization, thermalization and expansion timescales 
within irradiated samples that can be obtained with clus- 
ter experiments are needed for exploring the limits for ex- 
periments on single particle diffraction imaging [6-11]. 

In ref. [12] we investigated the ionization dynamics 
within Xe25oo clusters irradiated with rectangular vacuum 
ultraviolet (VUV) pulses of five different integrated pulse 



energy fluxes: F = 0.05,0.3,0.84,1.25,1.5 J/cm 2 . Pulse 
intensities were, I ~ 10 12 — 10 14 , W/cm 2 and pulse dura- 
tion was < 50 fs. For each flux the intensities and pulse 
durations were chosen in order to match the condition: 
I ■ At = F. Various intensities were tested at each fixed 
integrated radiation flux. The predictions obtained were 
then averaged over the number of events. This proce- 
dure enabled us to account for the non-linear response of 
the system to various pulse lengths and pulse intensities 
at higher radiation fluxes. This scheme followed the ex- 
perimental analysis: the experimental data were obtained 
after averaging the single shot data obtained with FEL 
pulses of a fixed integrated radiation flux. 

Here we study the non-equilibrium ionization dynamics 
within clusters of various sizes ranging from 20 to 90000 
atoms. These clusters are exposed to single rectangular 
VUV pulses of wavelength ^100 nm (photon energy ~ 
13 eV) and duration < 50 fs. Various pulse intensities 
are tested, however, the integrated radiation flux is kept 
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constant, F = I ■ At — 0.4 J/cm 2 . For comparison, if 
the pulse length is At — 50 fs, this integrated radiation 
flux can be achieved with the pulse intensity, I ~ 8 ■ 10 12 
W/cm 2 . The above pulse parameters correspond to those 
of the cluster experiment performed at the FLASH facility 
at DESY [13-15]. Results of this experiment show that 
ionization dynamics within the irradiated clusters strongly 
depends on cluster size. 

The present study applies the theoretical framework de- 
fined in ref. [12]. It completes the analysis performed in 
ref. [12] by evaluating the effect of the cluster size on the 
ionization dynamics. 

Our simulation tool solves the semiclassical kinetic 
equations for electron and ion densities within the irra- 
diated samples. Charges (represented as charge densities) 
interact with the mean electromagnetic field created by 
all charges and also with the laser field. Microscopic in- 
teractions enter these equations as rates. These rates are 
included into the two-body collision terms, and are esti- 
mated either from the experimental data or with theoret- 
ical models. Our model treats the following interactions: 
photoionizations, collisional ionizations, elastic scatterings 
of electrons on ions, inverse bremsstrahlung (IB) heating, 
shifts of energy levels within atomic potentials due to the 
plasma environment, and the shielded electron-electron in- 
teractions. 

This kinetic approach has two important features: (i) 
in contrast to hydrodynamic models, it can also follow the 
non-equilibrium sample evolution. Therefore it can deliver 
information on the timescales needed for reaching the local 
thermodynamic equilibrium (LTE) within a system, and 
(ii) it is computationally efficient for large samples. Fur- 
ther details of the method are discussed in refs. [12,16]. 
Based on the conclusions of ref. [12], we model the IB pro- 
cess with the enhanced IB rate and assume the plasma 
shifted energy levels of the atomic potentials. 

Below we discuss the results of our simulations per- 
formed for small and large clusters. The evolution of 
a sample hit by a single VUV pulse can be separated 
into two phases. The non-equilibrium ionization phase 
starts after the sample is exposed to the laser radiation 
and lasts until the saturation of ionizations from ground 
states is reached. Its duration may last up to several tens 
of femtoseconds depending on cluster size, pulse length 
(assumed to be < 50 fs) and pulse intensity (assumed to 
be ~ 10 12 — 10 14 W/cm 2 ). The semi-equilibrium expan- 
sion phase that follows after the ionization phase is much 
longer. Its timescale is of the order of picoseconds. Here 
we follow only the ionization phase and stop the simula- 
tion at entering the expansion regime. 

During the ionization phase electrons are released from 
atoms and ions. While they stay inside the sample, they 
are efficiently heated with the inverse bremsstrahlung pro- 
cess [17, 18]. Energetic electrons can then leave the cluster, 
creating an attractive Coulomb potential. Slower electrons 
are kept inside the sample, and the Coulomb attraction 
forces move these electrons towards the centre of the clus- 




Fig. 1: Schematic plots of charge distribution within an ir- 
radiated small cluster (left), and an irradiated large cluster 
(right) during the ionization phase. Positively charged outer 
shell coats the neutral core. Some electrons have escaped from 
the cluster (outer ionization). 

ter. Therefore, at the end of ionization phase the charge 
distribution within the sample shows a characteristic layer 
structure consisting of a neutral core and of a positively 
charged outer shell. This inhomogeneous charge structure 
was also described in refs. [12, 19-22]. In our case we ob- 
serve that the cluster core is dominated by ions of the high- 
est charges, however, the net charge of the core remains 
equal to zero. The core neutrality is due to the presence of 
quasi-free electrons bound within the core. The positively 
charged surface layer consists of ions of various charges, in 
particular, it contains ions of lower charges. 

We find that within smaller samples a large fraction of 
the electrons released during ionization processes is able 
to leave the sample early in the exposure. The width of 
the positively charged outer shell is then large, comparing 
to the radius of the neutral core (fig. Q] (left), figs. [2^- 
b). Therefore, we expect the Coulomb explosion to be the 
dominating mechanism of the overall cluster expansion. 

Within larger samples the number of energetic elec- 
trons that are able to leave the cluster during the ion- 
ization phase is small with respect to the total number of 
electrons released. Therefore, the width of the outer ion 
shell is small in comparison to the radius of the neutral 
core (fig. [T] (right), figs. [2j>d). It is expected that after 
the Coulomb-driven escape of the outer ions the core will 
slowly expand due to the hydrodynamic pressure of elec- 
trons and ions. Recombinations and ionizations (to and 
from excited states) will still occur within the neutral core. 
These processes may reduce the total ion charge within the 
core [12]. 

The Coulomb expansion is faster than the hydrody- 
namic expansion. The difference in these expansion rates 
can be useful at planning experiments on single particle 
imaging [9,10]. If a biomolecule could be coated with a 
layer of some other element or compound [10,23], only 
the outer part of the layer would undergo fast Coulomb 
explosion. The remaining neutral core containing the par- 
ticle under investigation would expand hydrodynamically, 
slower than the outer shell. Compared to the case of an 
uncoated sample, the coating would increase the radia- 
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Fig. 2: Formation of the outer shell of ions during the ion- 
ization phase. The charge density is defined as q ne t (r,t) = 
i-rii(r,t) — n e (r,t), where m (r, t) and n e (r, t) denote ion 
and electron densities. The maximal ion charge is, Nj = 8. 
Plots show the neutral core and the positively charged outer 
shell created within clusters of various sizes: a)-b) extensive 
outer shell for smaller clusters, c)-d) narrow outer shell for 
large clusters. Some electrons have escaped from the cluster 
(outer ionization) 

tion tolerance of the sample. This idea was investigated 
in ref. [10] for samples irradiated with X-rays. The sim- 
ulations performed within our study support the above 
scenario. 

Next we compare our results to the experimental data 
obtained with FLASH at DESY at 100 nm radiation wave- 
length [13,24]. The experimentally estimated histograms 
of ion charge and predictions on the average energy ab- 
sorption were obtained from the averaged time-of-flight 
(TOF) spectra. We remind here that the TOF detector 
could record charged particles (ions) only. There are no 
data on neutrals (atoms) available from these measure- 
ments. The experimental ion intensities were obtained by 
integrating and averaging the TOF signal over 100 subse- 
quent FEL pulses. Those intensities were not corrected for 
the relative geometric acceptances of the TOF detector or 
the MCP detector efficiencies for different charge states. 

First we present the simulation results for smaller clus- 
ters consisting on average of 20 and 70 Xe atoms. As 
stated previously, a large fraction of the electrons released 
during the ionization processes can leave the small clus- 
ters early in the exposure. The remaining electrons are 
not heated efficiently via IB processes due to their low 
density within the cluster. Only charge states up to +2 
(N — 20 atoms) and +4 (N = 70 atoms) are observed 
(fig. [3^i-b). The predictions obtained are in agreement 
with the experimental findings. The agreement is bet- 
ter for clusters of 20 atoms, for cluster of 70 atoms we 



Fig. 3: Ion fractions, Ni/N, within the irradiated xenon clus- 
ters at the end of ionization phase. The number, N, denotes 
here the number of Xe +l ions, and the number N is the to- 
tal number of ions. At t = fs these clusters contained: a) 
20 atoms, b) 70 atoms, c) 2500 atoms and d) 90000 atoms. 
For Xe20 cluster the experimental and theoretical results over- 
lap. For larger clusters b)-d) also the ion fractions within the 
positively charged outer shell are shown. The experimental re- 
sults are in agreement with the ion fractions from the positively 
charged outer shell. 

find a discrepancy between the predictions and the data. 
This discrepancy can be reduced, if we assume that the 
ions from the positively charged surface layer of the clus- 
ter give the dominant contribution to the experimentally 
measured ion charge spectra. This assumption has been 
discussed in detail in ref. [12]. The resulting outer shell 
predictions are shown in fig. [3p). Experimental data are 
within the model predictions. 

Further, we investigate large clusters built of 2500 and 
90000 xenon atoms. After the irradiation, only a small 
fraction of the electrons is able to escape from these large 
samples. The width of the outer shell created is small with 
respect to the radius of the neutral core. Following the 
discussion in ref. [12], in fig. [3p-d we plot the ion fractions 
obtained within the whole cluster and within its outer shell 
at the end of the ionization phase. Our predictions on the 
maximal ion charge are in agreement with experimental 
predictions, except of the absence of charges +8 within the 
largest cluster, N — 90000 atoms. The trend of the size 
dependence is correct: at fixed pulse energy the maximal 
ion charge created increases with the cluster size. 

We note that the ion fractions observed within the whole 
cluster at the end of the ionization phase will not corre- 
spond to those recorded by the detector during the exper- 
iment. Ions from the positively charged surface layer will 
be the first ones to escape from the sample, and they will 
reach the TOF detector with the unchanged charge state 
distribution. For the cluster core the situation is different. 
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Fig. 4: Average charge, Z, created within the whole irradiated 
cluster (red errorbars) and within the outer shell (green error- 
bars) as a function of the cluster radius, R. These estimates 
were obtained with pulses of different intensities and lengths 
but of fixed integrated radiation flux, F = 0.4 J/cm 2 , and 
then averaged over the number of pulses. Errorbars estimate 
maximal errors. Experimental data are plotted with stars. 



Fig. 5: Average energy absorbed per atom, E, within the ir- 
radiated cluster as a function of the cluster radius, R. Upper 
and lower limits (red and green errorbars) for the absorbed en- 
ergies are estimated within our model. These estimates were 
obtained with pulses of different intensities and lengths but of 
fixed integrated radiation flux, F = 0.4 J/cm 2 , and then aver- 
aged over the number of pulses. Errorbars estimate maximal 
errors. Experimental data are plotted with stars. 



At the end of ionization phase the core is a dense system 
of quasi-free electrons and ions, where recombination and 
ionization processes occur. During the long picosecond ex- 
pansion phase charges within the core have enough time 
for an efficient recombination. Therefore, the remnants of 
the core will be weakly charged or neutral. They will reach 
the detector at the end of the ionization phase, and then 
modify the ion charge state distribution recorded in the 
detector by increasing the participation of lower charges. 
This mechanism should be quantitatively verified with an 
expansion code, e. g. hydrodynamic code. This is, how- 
ever, beyond the scope of the present study. 

We focus now on the ion charge state distribution ob- 
served within the outer shell for large clusters, N = 
2500,90000 atoms. They are in good agreement with 
the experimental data (fig.[3b-d). This indicates that the 
recombination within the cluster core is efficient during 
the expansion phase. Also, the average charge predicted 
within the whole cluster and within the outer shell sup- 
ports this observation. This charge is plotted as a function 
of cluster radius in fig. The charge calculated within the 
outer shell follows the experimental results and starts to 
saturate for clusters of radius > 40 A, i. e. built of more 
than 3300 Xe atoms. 

Within small clusters, recombination processes will not 
be efficient due to the low electron density within the core. 
Therefore one can expect that the ion fractions recorded at 
the end of ionization phase within the whole cluster will 
approach those detected in experiment. Our results on 
small clusters of N = 20, 70 atoms support this scenario 
(fig.St,b). 

Below we show also the average energy absorbed per 
atom estimated with our model as a function of cluster 
radius (fig. [5]). With our model we can only obtain the 
upper and the lower limit for the absorbed energy. The 



upper limit assumes that during the further expansion of 
the sample no recombination processes are occurring. The 
lower limit gives an estimate of the energy absorbed per 
atom in case of the full neutralization of the sample during 
the expansion (full recombination). These limits are com- 
pared to the available experimental data on the average 
ejection energy per atom. Experimental data for larger 
clusters are within the model estimates. The energy ab- 
sorption per atom increases with the cluster size and sat- 
urates for larger clusters. This tendency is in qualitative 
agreement with the experimental data on the ejection en- 
ergies of three different ion charge states, +1, +4 and +7 
(fig. 3.11 in ref. [14]). 

We also list our predictions on the ionization and ther- 
malization timescales for electrons within the clusters. 
These timescales are measured from the start of the expo- 
sure at t = fs. We checked that if the cluster is irradi- 
ated with single rectangular VUV pulses of a flux fixed to 
I ■ At — 0.4 J/cm 2 , saturation of ionizations from ground 
states occurs on average after: tj ~ 2.5 — 3 -At. For longest 
pulses (At ~ 50 fs) tj ~ At and for shorter pulses (At ~ 4 
fs) tj ~ 4At. 

The thermalization timescale is defined as the timescale 
needed to achieve the LTE within a system. After this 
time the velocity distribution of electrons approaches the 
Maxwell-Boltzmann distribution. Within the large clus- 
ters only a small fraction of electrons can leave the sam- 
ple. The remaining electrons are then heated up in- 
side the cluster by the IB process, and they thermalize 
fast. As the non-equilibrium electron escape dominates 
the evolution of smaller clusters early in the exposure, the 
thermalization timescales for small clusters will be longer 
than for the large clusters. Thermalization timescales, 
tx, estimated at this particular integrated radiation flux, 
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F = I ■ At = 0.4, J/cm 2 are: i) for small clusters (20 and 
70 atoms) 0.8 • At < t T < 4 • At, and ii) for large clusters 
(2500 and 90000 atoms) 0.6- At < t T < 2.5- At, depending 
on the pulse length. For comparison, in ref. [12] the ther- 
malization timescale for the cluster Ve 2 5oo at the higher 
integrated radiation flux of F = 0.6 J/cm 2 was found to 
be > 0.3 • At. 

To sum up, we used a microscopic model based on 
the first principle Boltzmann approach to investigate non- 
equilibrium dynamics within atomic clusters of various 
sizes irradiated by single VUV pulses. This model is com- 
putationally efficient for small and large clusters, and in- 
cludes various processes that are relevant at VUV pho- 
ton energies. In particular, it uses the enhanced IB rate 
obtained using effective atomic potential as proposed in 
ref. [17, 18] and accounts also for the plasma environ- 
ment effects: the electron screening [25] and the lowering 
of the interatomic potential barriers (following the ideas 
of refs. [21,26,27]). Within our model cold electrons re- 
leased during photoionizations are heated by the IB pro- 
cess. Highly charged ions are created during collisions of 
energetic quasi-free electrons with ions. Electrostatic in- 
teractions among charges, interactions with the external 
laser field and shielded electron-electron interactions are 
treated. Recombination processes are neglected due to 
their low rates expected and short simulation timescales 
(< 100 fs) [12]. 

The results obtained demonstrate the different ioniza- 
tion dynamics of small and large clusters. Predictions ob- 
tained at various cluster sizes on: (i) the maximal and av- 
erage ion charge, (ii) ion fractions within the outer shell, 
and (iii) average energy absorbed per atom were compared 
to the experimental data. They were in good agreement 
with these data. However, a quantitative verification of 
the outer shell predictions should be performed with an 
expansion code. 

We have shown that our non-equilibrium Boltzmann 
solver is a useful tool for studying the evolution of irra- 
diated samples. Several interesting problems may be in- 
vestigated in future with this code and with an expansion 
code, e. g. the evolution of clusters irradiated with XUV 
and X-rays, the evolution of mixed clusters or mechanisms 
of slowing down the cluster explosion. 

Also, if ever ultrashort FEL pulses are available (1 < 
At < 10 fs) [28-31], the evolution of the samples irradi- 
ated with these pulses may be fully non-equilibrium during 
the exposure. As the applicability of the semi-equilibrium 
hydrodynamic models at these pulse lengths is much re- 
stricted, we expect that the non-equilibrium Boltzmann 
solver can be particularly useful for studying the prop- 
erties of larger samples irradiated with ultrashort FEL 
pulses. 
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